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Abstract 

We report on the calibration of a superfluid 3 He bolometer developed for the search 
of non-baryonic Dark Matter. Precise thermometry is achieved by the direct mea- 
surement of thermal excitations using Vibrating Wire Resonators (VWRs). The 
heating pulses for calibration were produced by the direct quantum process of quasi- 
particle generation by other VWRs present. The bolometric calibration factor is 
analyzed as a function of temperature and excitation level of the sensing VWR. 
The calibration is compared to bolometric measurements of the nuclear neutron 
capture reaction and heat depositions by cosmic muons and low energy electrons. 
The comparison allows a quantitative estimation of the ultra-violet scintillation rate 
of irradiated helium, demonstrating the possibility of efficient electron recoil event 
rejection. 
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1 Introduction 



Superfluid 3 He-B at ultra-low temperatures is an appealing target material 
for bolometric particle detection [T|f2|l3"] . with complimentary features to the 
currently most performing germanium- and silicon-based detectors like Edel- 
weiss [4] and CDMS [5J. Among the attractive features of 3 He are the clear 
signature of neutron background events due to the nuclear capture reaction 
of neutrons in 3 He. Furthermore, the huge density of unpaired neutrons leads 
to a significantly enhanced axial interaction cross-section within a large class 
of WIMP models [6f7] . Its single thermal energy reservoir with the possibility 
of direct thermometry, the very low experimental base temperature of about 
100 /iK, together with the possibility of electron recoil rejection make 3 He a 
very promising target material for Dark Matter search. 

The nuclear neutron capture reaction 3 He+n — > 3 H+p releases 764 keV kinetic 
energy to its products, and was detected bolometrically in the superfluid [3]. 
The comparison of the detected neutron peak to bolometric calibrations of 
the detector was then used [S] as a test of topological defect creation scenarios 
in fast second order phase transitions [H] . A bolometric calibration method of 
the detector was first described in [10]. In the following years, the detection 
threshold of the bolometer could be lowered, which allowed the identification 
of a broad peak at about 60 keV, that could be attributed to cosmic muons, in 
agreement with numerical predictions. In very recent measurements on a Dark 
Matter detector prototype, the detection threshold reached the keV level and 
the low energy electron emission spectrum from a 57 Co source could be resolved 
|ll|fT2] . Since a neutralino interaction with a 3 He nucleus in the detector is 
expected to deposit a maximum energy of about 5.6 keV, the current detector 
design therefore already displays the required bolometric sensitivity. 

The recent improvement of the sensitivity was partly made possible by low- 
ering the working temperature of the detector cells from about 150 to 130 
/iK. This temperature decrease of less than 20% represents a decrease in the 
thermal enthalpy of the cells by an order of magnitude. On the other hand, 
at lower temperature the weaker coupling of the VWR to the superfluid re- 
sults nevertheless in a greater response time of the thermometer. Parallely, 
the non-linear velocity dependence of the friction with the superfluid [13], is 
also of greater importance at the lower temperatures and the higher VWR- 
response signals currently used. This article hence proposes a generalization 
of the methods described in [TU] for a more profound and comprehensive un- 
derstanding of the bolometric detector (sections 2 to 4) and its calibration 
by mechanical heater pulses (section 5). The comparison of the bolometric 
calibration with detection spectra from known energy sources like neutrons 
(section 6) as well as low energy electrons and muons (section 7) allows then 
a precise estimation of the ultra-violet (UV) scintillation rate of helium for 
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several types of ionizing radiation. 



2 Principle of detection 

The current 3 He particle detector [TT][T2"] consists of [TO] three adjacent cylin- 
drical copper cells, of volume V = 0.13 cm 3 each, immersed inside a superffuid 
3 He bath at ultra-low temperature, at about 130 /xK. A 200 diameter ori- 
fice in the wall of the cells connects the 3 He inside the cell with the surrounding 
heat reservoir of superffuid. An energy deposition by an incident particle cre- 
ates a cloud of quasiparticles in the cell, which comes to internal thermal 
equilibrium via collisions on the cell walls within ~ 1 ms. The excess ballistic 
quasiparticles, of momentum p w pp, then leak out of the cell with a time 
constant t& of a few seconds, which is determined mainly by the size of the 
orifice . 



Fig. 1. Schematic setup of the 3-cell bolometric detector prototype. Each cell con- 
tains a VWR thermometer. The cells are in weak thermal contact with the outer 
3 He bath through the orifice. The presence of three adjacent cells allows discrimina- 
tion through coincident detection. One of the cells contains the Co source (section 
7), another one contains an extra VWR for the purpose of bolometric calibration 
by heater pulses (section 5). 

The Vibrating Wire Resonators (VWRs) used for thermometry are thin (4.5 
/im diameter in our setup) superconducting NbTi filaments bent into an ap- 
proximately semi-circular shape of a few mm, with both ends firmly fixed |14j . 
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The VWR is driven by a Laplace force imposed by an AC current close to its 
mechanical resonance frequency (~500 Hz), and oscillates perpendicularly to 
its main plane with an rms velocity v. The motion is damped by frictional 
forces, of total amplitude F(v), mainly due to momentum transfer to the 
quasiparticles of the surrounding superffuid. Relatively fast thermometry is 
then achieved by measuring continuously the damping coefficient W ex F/v 
of a VWR driven at resonance. In the low oscillation velocity limit, the friction 
is of a viscous type F oc v . W is thus velocity independent and can be written 
as [15] 

Wo = a exp{-A/k B T), (1) 

where A is the superffuid gap at zero temperature. The value of the prefactor 
a is determined by both the microscopic properties of the liquid and the 
geometry of the VWR. At bar and for the VWRs used in the detector, of 
density p = 6.0 g/cm 3 and radius a = 2.3 /xm, a is of the order of 1-2 10 5 Hz, 
depending on the exact geometrical features. In quasi-static conditions, the 
low velocity damping coefficient Wo is measured as the Full Width at Half- 
Height A/2 of the mechanical oscillator's lorentzian resonance in frequency 
space. 

At finite oscillation velocities however, effects beyond the two-fluid model dis- 
covered in Lancaster, lead to a velocity dependent friction which is actually 
decreasing at higher velocities (Fig. [2]). A model proposed by Fisher et al. 
[T3~]n5] leads to a velocity dependent damping 

W L (v) = W x L(v/v ), (2) 

where the function 

L{u) = (1 - e~ u )/u (3) 

gives a very good approximation of the reduction of the damping of the VWR 
for v up to a characteristic oscillation velocity Vq ~ ksT/pp « 1-2 mm/s. 

At higher velocities, of about 3-4 mm/s, the local flow field around sharp edges 
on the VWR surface starts to get turbulent [16J. At even higher velocities, 
the critical velocity for pair-breaking is achieved. Both cases correspond to 
a strong enhancement of the friction, with a strongly non-linear, sometimes 
even discontinuous, velocity dependence (Fig. [2]). While the Lancaster- type 
non-linear coupling to the superfluid can easily be accounted for by using 
equation (j2J), the velocity has thus still to be kept small enough for both the 
locally turbulent and the pair-breaking regime to be avoided. 
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Fig. 2. Isothermal force- velocity (= (I,V)) curve of a VWR immersed in super- 
fluid 3 He at ultra-low temperatures. For rms oscillation velocities below 2 mm/s, 
the damping is well described by ([2]). The decrease of the damping coefficient (~ 
inverse slope) at higher velocities actually leads to an increase of the calibration 
factor (section 4). At velocities above ~2.5 mm/s, the damping increases again. 
Points 1 and 2 show two discontinuities of the force- velocity curve (see text). The 
pair-breaking velocity is reached at about 4.5 mm/s. 

3 Response of the detector to a heat release 



As the quasiparticle density, the temperature, and thus W as denned by 
([1]), increase nearly instantaneously (~ 1 ms) after a particle impact, the 
mechanical VWR only responds to this modification over a timescale inversely 
proportional to the dissipative coupling 

1 

T~w 777" 
7lW L 

This can therefore result in non- negligible response times (>ls) for low damp- 
ing coefficients, i.e. narrow resonances and low temperatures. While recording 
transitory heating events, the dynamically measured damping W mes (t) is thus 
a non-equilibrium measurement of W(t). 

A simple model of a sudden increase of W from its baseline W\, ase by an 
amount A <C Wi, ase at t = 0, and the subsequent exponential relaxation of the 
quasiparticle excess via the orifice 



(4) 



W{t) = W base + Ae- t/Tb Q{t), 



(5) 



leads to a dynamic damping measurement given by 



W mes {t) = W base + A- 



n 



n 



©(*)■ 



(6) 



G(t) is here the Heaviside function accounting for the instantaneous heat input 
and since A is a small perturbation, t w « constant = l/itWbase is assumed. 
The good agreement of (El) with experimentally detected events can be seen 
on figure El The maximum geometrical height H = \W mes — Wb aS e\max of the 
response peak is then related to A by the function 



G(W base ,T b ) = H/A=(T w /r b )-> 



b-T W 



(7) 



For say W base = 0.414 Hz (T = 129 /jK) and r b = 5.0 s, the slowing down 
of the VWR response is hence responsible for the loss of 29 % of the signal 
amplitude. 
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Fig. 3. Bolometrically recorded event (~ 100 keV) (circles). The fit by © uses 
the input parameters W base = 414 mHz and thus r w = l/irW base =0.77 s. The cell 
relaxation time r b = 5.0 s is fixed for all events of the entire experimental run. The 
only free parameters are the start time of the event and and A. Here the fit fields 
A = 56.0 mHz, in good consistency with H x G~ 1 (W base , r b ) = 55.4 mHz. 
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4 Calibration factor 



Since both the enthalpy and Wo are roughly proportional to the quasiparticle 
density, the bolometric calibration factor a of a given cell, defined as the 
conversion factor relating the amplitude A (<C Wo) of the rise in Wq in (jSJ) to 
the heat U released 



a = A/U, (8) 

is a slow function of temperature. Around 130 fiK, the temperature depen- 
dence of the specific heat of 3 He-B is given to a good approximation by [Toll?] 



C = Co (T c /Tf 2 exp(-A/k B T), (9) 

where T c is the superfluid transition temperature and Co = 1.7 mJ/Kcm 3 . 
Writing 5U = CV5T allows to express the calibration factor using (JTJ), (jHJ) 
and <M) in the low velocity limit as 



cto(T) = -L (10) 

k B CoVT? 12 VT 

Taking into account the velocity dependence of the damping (j2]) leads to an 
enhancement of the calibration factor at higher velocities over <Jo by a factor 



/(«,«') = <7(„)/<7„ = «' 



sinh(w/2) 



e 



-u/2 



sinh(^2) 



where t' = ksT/A and w = v/v . At velocities f ~ v , this results consequently 
in an enhancement of the calibration factor by about 12 %. This analysis is 
valid for velocities up to ~ 1.5 vq, before the other dissipative mechanisms 
discussed set on. 

Adding in (JJJ) leads thus to the conclusion that the relation of the geometric 
height H of a peak following a heat deposition U depends on both the baseline 
Wbase and the rms velocity v of the oscillation following 



H = a f(u, t') x G(W base , n) x U. (12) 
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5 Bolometric calibration by heater pulses 



At ultra-low temperatures, simulation of heating events by Joule heating is 
inefficient in superfluid 3 He because of the diverging thermal (Kapitza) re- 
sistance at the solid-liquid boundaries. Bradley et al. [3] proposed a heating 
method of the superfluid, based on the mechanical friction of the oscillating 
VWR with the liquid. Energy is injected through a second VWR present (the 
" heater" VWR, as in the lowest of the 3 cells in Fig. [Q , by driving it over a 
short time with an excitation current at its resonance frequency. The energy 
is then dissipated to the superfluid via the velocity dependent frictional cou- 
pling. One should bare in mind that the electrically injected power is firstly 
transformed into mechanical energy of the resonator before being transfered 
to the liquid by friction. Therefore, even for a short excitation pulse of the 
heater (St ~ 100 ms), the heat release to the liquid nevertheless takes place 
on a characteristic timescale given by (j4]), where the damping W is relative to 
the heater wire. 

The energy deposition by a heater pulse is hence not instantaneous, even in 
the case of extremely short pulses. We have estimated quantitatively this heat 
release, in order to verify wether a mechanical heater pulse is indeed an ac- 
ceptable 'simulation' of a real - instantaneous - event. If one assumes velocity 
independent damping of the heater resonator, its amplitude of oscillation must 
increase linearly during the pulse and then relax exponentially with timescale 
t. We have modeled numerically under these assumptions the heat release to 
the superfluid in such a heater pulse, as well as the response of the thermome- 
ter wire in the bolometric cell. One result and its comparison to the response 
calculated for an instantaneous energy release as say by an incident particle, 
are given in Fig. |H One can see that even though the shape of the thermome- 
ter's response is quite different, especially during the rise of the response, the 
maximum amplitude is still the same within less than 1 % for Wb ase > 0.3 Hz. 
A fit of the response by ([6]) is therefore not appropriate in the case of mechan- 
ical heater pulses, while the geometrical maximum height H of the response 
peak is still a good parameter for comparing pulses to instantaneous events. 

A more serious issue for the bolometric calibration by mechanical heater pulses 
are the intrinsic losses within the heater wire. The intrinsic damping (i.e. in 
the absence of liquid) of a VWR can be well modelled by a viscous friction 
term which therefore simply adds as a constant extra term Wi nt to the mea- 
sured damping W in ([I]). In our experiment, we measure the heater VWR 
intrinsic width as the limiting value of the damping in the low temperature 
limit, to be W^ t ater = 77 ± 6 mHz. As the thermal quasiparticle damping 
vanishes at low temperatures, the intrinsic losses in the heater represent a 
fraction W^ t ater / (W^ t ater + W^ eater ) of the injected energy and constitute a 
non-negligible correction to the energy injection through the heater wire. 
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Fig. 4. Simulation of the response of the thermometer VWR to an instantaneous 
heat deposition at t = 1 s (dotted line) and to a mechanical heater pulse of same 
total energy but releasing the heat slowly (continuous line) because of T^, eater . While 
the response to the instantaneous event is well described by ([6]), the response to the 
slow heat release has a different shape during the rise, but is then rather identical 
exept for a time shift due to the slower release of the energy by the heater. The 
height H of the two peaks is nevertheless identical within 1 % at 0.3 Hz (128 fiK) 
and much less at higher temperatures. 

Let us now consider the energy and velocity dependence of the detector re- 
sponse. At constant oscillation velocity, the height H of the response is exper- 
imentally observed to be linear in the injected energy in the range of a few 
hundred keV (Fig. |3j), which confirms the assumptions leading to (jHJ). Fur- 
thermore, a comparison of the observed slopes versus the VWR's oscillation 
velocity (Fig. [6]) are in good agreement with the enhancement of the response 
at higher velocities, as given by (1111) . This analysis allows therefore to extract 
the velocity independent limit of the calibration factor a . 

Taking into account intrinsic losses in the heater as well as the features of the 
thermometer VWR response deduced in the previous sections, one can fit the 
mechanical pulse calibration data versus baseline (i.e. versus temperature) as 
shown in Fig. If W^ t ater is left as a free parameter, the fit yields a value 
of W^ t ater = 74 mHz and a calibration factor of cr = 3.8 10~ 4 Hz/keV at a 
baseline of 1 Hz (JTj). Imposing Wj^ t ater = 77 mHz leaves the overall prefactor 
of Co (equivalent to the knowledge of a) as the only free fitting parameter, 
with a value only 0.5% different. 
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Fig. 5. Height H of the response of the thermometer VWR to mechanical heater 
pulses of energy E (bullets) and linear fit (data taken at Wbase = 1-7 Hz). 
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Fig. 6. Dependence of the observed slopes of H(E) as in Fig. [5] as a function of 
the reduced oscillation velocity of the thermometer VWR {Wbase = 4 Hz). The 
continuous line represents the expected evolution from (llip , the only free parameter 
being H/E(v = 0). 
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Fig. 7. Calibration of the observed H/E in the zero velocity limit from mechanical 
heater pulses (•). The data are then corrected by ([7j) to account for the slowing down 
of the thermometer at low temperatures (o). Considering the fact that a fraction 
Wi£t' ter /Wto t ater of the injected energy is dissipated inside the heater VWR, the 
result can be fitted by (|10|) (dashed line). The only free parameter is the overal 
prefactor of (|10p which allows then to express &o(T) (continuous line). 
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6 Neutrons 



The heat release to the superfluid following a nuclear neutron capture reaction 
3 He + n — > 3 H + p was observed in [8] to be at bar about 15 % less than the 
764 keV known to be released by the reaction. For a microscopic calculation 
of the stopping of the neutron capture products in 3 He see [18] . A substantial 
part of the energy deficit can be attributed to ultra-violet (UV) scintillation 
[TU] . The UV fraction was estimated in [5] to account for about half of the 
energy deficit, concluding that the remaining ^60 keV were trapped in the 
form of topological defects. 

Neutrons are an extremely valuable tool for the characterization of the super- 
fluid bolometers because they provide a narrow and cell-geometry independent 
peak at rather high energies. Once the bolometric heat release of neutrons in 
3 He inambiguously determined, no other calibration, e.g. by pulses, is nec- 
essary. The comparison of the evolution of the observed neutron peak with 
temperature also allows to test directly, without the use of mechanical heater 
pulses, the assumptions leading to f[T2"|) . 
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Fig. 8. Detected events in the energy range of a few hunderd keV in presence of the 
AmBe neutron source, as a function of Wb ase . Note that several points do not ac- 
tually represent single neutron capture events, but the sites of larger concentration 
of events at a given Wb ase (circles) correspond to the neutron peak at that temper- 
ature. The continuous line represents the temperature dependence of the expected 
height H corresponding to the neutron peak as given by (|12l) . 

The dependence of the neutron peaks' position on the baseline (i.e. on tem- 
perature) is shown on Fig. El The agreement of the baseline dependence as 
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given by (TT2T) is good except for the very lowest temperatures where an ab- 
solute deviation from (fl2|) of about 6 % is observed at Wb ase =100 mHz. The 
analysis yields a bolometric heat deposition of 652±20 keV at bar, in very 
good agreement with earlier measurements, where 650 keV were found under 
the same conditions [8]. Wether the observed energy deficit of about 110 keV 
can be significantly accounted for by topological defects is nevertheless still an 
open question. An extensive study of the pressure dependence of this deficit 
is in progress and might give some new insight on this question [20J. 



7 Scintillation yield in helium of muons and low energy electrons 

The detection of electrons of 7 and 14 keV emitted by a 57 Co source within one 
superfluid 3 He-B bolometer has been reported recently [12]. The comparison 
of the detected electron spectrum with the bolometric calibration factor of 
that cell, determined as described in sections 2-5, yields a deficit of 25±5 % 
in the bolometrically detected energy. We therefore define the scintillation 
yield in the following as the fraction of the injected energy that is released as 
radiation. 

In helium, the distance between two consecutive collisions of electrons in the 
keV range is a fraction of a /im. This distance is much larger than that in the 
case of a-particles or the neutron capture products. Upon electron irradiation, 
the energy is hence deposited with much lower density and at no time the 
density of heat is such that the superfluid could heat up to the normal state. 
In the case of electrons, the energy deficit is therefore entirely due to UV 
scintillation. The scintillation yield is very sensitive to the density of the energy 
deposition. Our results can be compared to measurements by Adams et al. [2T] 
who find a 35 % UV scintillation yield in 4 He upon irradiation by electrons 
of a few hundred keV. Nevertheless, as later emphasized by McKinsey et al. 
[i~9f22] . this measurement only covered the fast emission pulse of the first 10- 
20 nanoseconds. A substantial part of the UV emission, estimated to ~ 50 % 
of the fast pulse contribution by these authors, takes place at much later 
times, which brings the UV scintillation yield to a total of about 50 % for high 
energy electrons. At lower energies on the keV level, this fraction is expected 
to decrease rapidly. Our findings of the scintillation rate at ~ 10 keV represent 
thus a very reasonable low energy extension of the measurements by McKinsey 
et al. 

The muon events also display in our detector a clearly resolved spectrum. 
Figure [9] shows data taken from [23J, comparing numerical simulations to 
measurements taken within the MACHe3 collaboration on the muon peak. 
A GEANT 4.0 code computer simulation of the muon energy release in our 
3 He cell, integrated over a large spectrum of incident muon energies, leads to 
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a broad peak at about 67 keV. The dispersion of the energy distribution is a 
result of the muons tracks crossing the cell geometry in all space directions 
with the known angular dependence of radiation intensity [23] . The muon en- 
ergy release to the target material is rather independent on the energy of the 
incident muons in the range from 1 to 10 GeV covering the maximum of the 
energy distribution of muons at ground level. 

The width of the experimentally detected muon spectrum corresponds to the 
intrinsic resolution of our detector (3 % ) combined to the geometrical broad- 
ening. The muons produce exclusively electron recoils in the cell producing 
delta rays, i.e. electrons scattered from the atomic shells, energy-distributed 
following E~ 2 . For one incident muon, roughly 10 of such electrons have an 
energy greater than 250 eV and only 5 greater than 1 keV. 




In our experiment we have collected data in the muon energy range (Fig. [9]) 
during 19 hours. A broad maximum of events with an average energy of 52 
keV is clearly seen. The comparison between the bolometric measurement and 
the simulation of the energy input by muons therfore yields again an energy 
deficit of about 25 %. The missing energy being the result of the scintillation 
of UV photons as for incident electrons, this observed deficit again allows to 
quantify the scintillation rate of helium for this type of irradiation. 
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Table 1 

Measured energy release by different particles and their scintillation rates 



particle 


muons 


electrons 


p+ 3 R 


initial kinetic energy (keV) 


2-4 10 6 


14 


571+193 


energy release to target 3 He (keV) 


~67 


14 


764 


energy detected bolometrically (keV) 


~52 


11 


652 


scintillated fraction 


0.25 


0.25 


< 0.15 



8 Conclusions 



We have presented a detailed description of the method of bolometric calibra- 
tion of the detector cells based on superfluid 3 He, by mechanical pulses. This 
description provides a general and comprehensive picture of the bolometer 
under different conditions of use (excitation level, temperature). 

The results of the calibration are compared to measured heat depositions by 
the nuclear neutron capture reactions as well as muon impacts and low energy 
electron irradiation. A deficit of about 15 % is found in the case of neutrons, 
in good agreement with previous measurements at bar j8] , in which case this 
deficit is partly associated to vortex creation. In the case of high energy muons, 
as well as electrons in the lOkeV range, a deficit of about 25 % is found which 
can be entirely attributed to UV scintillation emission. It is not surprising to 
find the scintillation rates resulting from these two types of irradiation to be 
of the same order since the much larger incident energy of cosmic muons is 
compensated by their larger mass. 

Now that our detector has achieved the required sentitivity for resolving recoils 
in the keV energy range, we focused on the the feasibility of a electron recoil 
rejection by parallel scintillation detection. Since we found the UV scintillation 
yield not to be small for both electrons and muons down to the energy range 
of interest, this result gives experimental evidence that the parallel use of 
a scintillation detector would allow to reject efficiently electron-, muon- and 
most likely 7-contamination in a bolometric detector based on superfluid 3 He 
for the search of non-baryonic Dark Matter. The optimum design of such a 
parallel scintillation detector will be discussed in a future work. 
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